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ABSTRACT 



Context. White dwarf spectra have been widely used as a calibration source for X-ray and EUV instruments. The in-flight effective area 
calibration of the reflection grating spectrometers (RGS) of XMM-Newton depend upon the availability of reliable calibration sources. 
Aims. We investigate how well these white dwarf spectra can be used as standard candles at the lowest X-ray energies in order to gauge the 
absolute effective area scale of X-ray instruments. 

Methods. We calculate a grid of model atmospheres for Sirius B and HZ 43A, and adjust the parameters using several constraints until the 
ratio of the spectra of both stars agrees with the ratio as observed by the low energy transmission grating spectrometer (LETGS) of Chandra. 
This ratio is independent of any errors in the effective area of the LETGS. 

Results. We find that we can constrain the absolute X-ray spectrum of both stars with better than 5 % accuracy. The best-fit model for both 
stars is close to a pure hydrogen atmosphere, and we put tight limits to the amount of helium or the thickness of a hydrogen layer in both stars. 
Our upper limit to the helium abundance in Sirius B is 4 times below the previous detection based on EUVE data. We also find that our results 
are sensitive to the adopted cut-off in the Lyman pseudo-continuum opacity in Sirius B. We get best agreement with a long wavelength cut-off. 
Conclusions. White dwarf model atmospheres can be used to derive the effective area of X-ray spectrometers in the lowest energy band. An 
accuracy of 3-4 % in the absolute effective area can be achieved. 



^ ^ Key words. Stars: atmospheres - white dwarfs -X-rays: stars - X-rays: general 



1. Introduction 

DA white dwarfs such as Sirius B and HZ 43A can be used 
very well as calibration sources for UV and X-ray spectrom- 
eters. This is because their spectra are simple and dominated 
by hydrogen features. In some cases there are small traces of 
He visible in the spectrum. The simplicity of their constitution 
and in particular for the hotter stars the straightforward physics 
of the atmosphere allows ab initio calculations of the emerging 
spectra. The models depend only upon a few free parameters, 
predominantly the effective temperature T^ff, the surface grav- 
ity g and the relative helium to hydrogen number density. The 
effective temperature and surface gravity can be deduced from 
optical or UV spectra (Balmer and Lyman series). Despite the 
relatively low surface temperatures, that are generally below 
10^ K, DA white dwarfs emit soft X-rays because of the low 
opacity of the atmosphere for the high energy photons that are 
produ ced in the deep and hot inner layers of the star (Shipman, 
1976h . Therefore, if an independent estimate of the radius of the 



star as well as its distance are known, the ab initio calculations 
yield in principle a model spectrum for the full spectral range 
including the X-ray band, with known shape and normalisation 
(excluding ISM opacity). By comparing such a model spectrum 
to the X-ray spectrum as observed with a given instrument, the 
effective area of that instrument can be calibrated. 

Even when the distance or radius of the star are not accu- 
rately known, measurements of the absolute optical flux of the 
star can be used to scale the X-ray flux of the star In the op- 
tical band, absolute fluxe s can be determined with an accuracy 
of the order of a percent ( Holberg & Bergeronl2006t) . 

The Chandra low energy transmission grating spectrometer 
(LETGS) is currently the most sensitive high-resolution X-ray 
spectrometer operating in the softest X-ray band. The effec- 
tive area calibrati on as produ ced by the Chandra X-ray Center 
was described bv lPease et al.l ( 12000,) : essentially, a pure hydro- 



gen, non-LTE model for Sirius B with T^ff = 25 000 K and 



g = 10 ms was used dHolberg et al.lll998r) . with an esti 
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mated flux uncertainty of less than 10 %. A correction fac- 
tor to the effective area was determined by comparing the ob- 
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served spectrum of Sirius B to this model. Using the result- 
ing corrected effective area, the observed LETGS spectrum of 
HZ 43A agreed within 10-15 % with a model s pectrum for 



that source. Later refinements of the effective area dPease et al 



2003h mainly concerned the mid-energy range, based on spec- 
tra of the blazars PKS 2155-304, 3C 273 and the isolated neu- 
tron starRXJ 1856.5-3754. 

Alternatively, at SRON, J. Kaastra and J. Heise derived the 
effective area of the LETGS following a different approach de- 
scribed in an internal report (Kaastra 2000). Basically, a grid 
of models was calculated us ing a version of Tlusty dHubenvl 
19881: iHubenv & Lanzl[l995 l) available at that time, with the 
effective temperature and gravity as free parameter, for both 
Sirius B and HZ 43A. A 6-dimensional grid search was done 
(varying for each source the effective temperature, gravity and 
interstellar absorption column), in order to find the best match- 
ing spectrum using as constraints (i) that the ratio of both model 
spectra should match as closely as possible the measured count 
ratio with the LETGS; (ii) that the 70-170 A spectrum of 
Sir ius B must be within the error bars of the model as derived 
by iHolberg et al. (11998 !): and (in) that the model parameters 
should not be too far off from the parameters for both stars as 
listed in the literature. 



Th e ef fective area of the LETGS as derived by iKaastra 
( I2OOOI) and lPease et aD ( I2OOOI) differ typically by 10 %, but the 
differences are not constant as a func ti on of wavelength. 

More recently, Beuermann et alj (I2OO6 ) have tackled the 
problem of cross-calibration again using LETGS spectra of 
Sirius B, HZ 43A as well as RXJ 1856.5-3754. Assuming a 
double blackbody model for the latter source, and running grids 
of white dwarf models for the first two sources, they derived 
the effective area of the LETGS by fitting simultaneously the 
parameters of these stars and the effective area correction fac- 
tor. Based on their solution, they concluded that the Extreme 
Ultraviolet Explorer (EUVE) short wavelength (SW) detector 
effective area is off by ~ 15 % and that the ROSAT PSPC de- 
tector agrees within a few percent with the LETGS. 

However, the assumption of a double blackbody model for 
the neutron star is not unquestionable. For instance, the source 
could have a more complicated atmosphere, or there may be 
multiple hot spots or a single spot with multi-temperature struc- 
ture. Als o, comparing our m odel calculations to the models de- 
rived by iBeuermann et al. 1 12OO6) shows some significant dif- 
ferences, even if the same stellar parameters are used (see 
Sect. 15.5b . This gives us sufficient motivation to re-examine 
carefully the calibration of the LETGS and the model spectra 
employed for both white dwarfs. 

This paper is the second of a series of three intended to de- 
rive the absolute effective area of the reflection grating spec- 
trometer (RGS) of XMM-Newt on and through this of other 
instruments. In our first paper (Kaastra et al. 2008ah we have 
studied the RGS spectrum of Crab and derived accurately the 
interstellar absorption towards that source. However, the main 
uncertainty of 10 percent on the absolute flux of Crab could 
not be resolved. In the present paper we show how we can ac- 
curately calibrate the low-energy effective area o f the Chandra 



results using blazar spectra taken simultaneously with the RGS 
and LETGS. This leads to a reduction of the uncertainty in the 
Crab flux and absolute effective area of the RGS to about 3 
percent. 

2. Data analysis and spectral modelling 

We follow the same procedures as iKaastral (I2OO0I) but use the 
currently best available spectral models and boundary condi- 
tions. In Sect. 12. II we derive the intensity ratio of the spectra of 
Sirius B and HZ 43 A as measured by the Chandra LETGS. 
We will adjust the parameters of both stars using a number 
of boundary conditions until the best agreement with this ob- 
served intensity ratio is obtained. Sect. l2.2l describes our model 
for the white dwarf atmospheres. These model spectra essen- 
tially give the surface flux of the atmosphere, so we describe 
the scaling procedure to obtain the flux received at Earth in 
Sect. 12.31 In Sect. 12.41 we list all constraints that we apply to 
our models, and Sect. 12.51 describes our spectral fitting proce- 
dure. 

2.1. Data analysis 

The data used in this paper are summarised in Table [T] All 
data were obtained with the high resolution camera (HRC-S) in 
combination with the low energy transmission grating (LETG). 
We will call this combination here LETGS, as in our paper we 
do not use the LETG combined with the ACIS-S detector. 

For Sirius B we used all three LETGS spectra taken 
October 26-29, 1999. The net total exposure time is 44.6 ks. 
Recently (18 January 2008) the source was observed again for 
calibration purposes but as the source was close to the detector 
edge we do not use these data. For HZ 43A we used observa- 
tion ID 59, observed November 12, 1999, with a net exposure 
time of 38.9 ks. HZ 43A is monitored regularly since then, but 
as we are only interested in the ratio of the spectra of Sirius B 
to HZ 43 A, and HZ 43 A is the strongest source, the statistical 
uncertainty on the ratio is dominated by the statistical uncer- 
tainty of the Sirius B spectrum. More over, the sensitivity of th e 
LETGS slowly decreases over time jBeuermann et al. 20061) . 
and therefore it is important to compare observations not too 
far apart in time. 



The data processing is similar to iKaastra et al.l (120021) . 



LETGS using white dwarfs. In our third paper ( IKaastra et al 



2008bh we combine these results on white dwarfs with the Crab 



Periods with high background are filtered out in two steps. In 
the first step, the observation is split into 1 s intervals; all in- 
tervals where the total detector count rate saturates (more than 
180 counts s"') are discarded. After this a light curve for the 
zeroth order spectrum of the source is created, binned in 100 s 
intervals. This serves to check the constancy of the source. All 
intervals with less than 50 s exposure are discarded. The re- 
maining exposure time is then corrected for dead time using 
the standard dead time correction factors taken every 2.05 s, 
provided with each dataset. 

A few other health checks are made; in one of these, we 
compare the average pulseheight (PHA) of the dispersed spec- 
trum with a reference spectrum. This comparison is done in 
2 - 20 A wide bins, and the average PHA difference of the 
spectrum compared to the standard is calculated (see Table[T]i. 
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Source 


Sirius B 


Sirius B 


Sirius B 


HZ43A 


ObsID 


1421 


1452 


1459 


59 


start date 


28 Oct 1999 


26 Oct 1999 


27 Oct 1999 


12 Nov 1999 


Duration (ks) 


25.3 


28.0 


12.0 


40.2 


Net exposure (ks)" 


13.8 


19.0 


11.8 


38.9 


0th order counts/s 


1.282 ±0.010 


1.249 ±0.008 


1.310±0.011 


8.826 ±0.015 


PHA difference* 


-5.6 ±2.0 


-16.2 ± 1.9 


-6.5 ± 2.0 


-7.0 ± 1.4 



a Dead-time corrected exposure time after filtering out periods of high background. 

b The PHA (pulseheight) difference is the average PHA of the spectrum relative to a standard spectrum. 



Clearly, the average PHA in Sirius B observation 1452 is sig- 
nificantly lower than for all other observations. Also, the av- 
erage zeroth order count rate is lower: while the count rates 
for observation 1421 and 1459 are consistent with an average 
value of 1.295 ±0.007 counts s"', for observation 1452 we have 
1.249 ± 0.008 counts s"', or 3.6 % less counts. A smaller effect 
is present in the dispersed spectra of Sirius B. We have fitted 
all three spectra individually and found that their shapes are - 
within the error bars - consistent with the same model, but the 
average flux for observation 1452 is 0.8±0.4 % smaller than for 
the other observations. As observation 1452 contributes 43 % 
to the total exposure time for Sirius B, we kept the data but cor- 
rected the fluxes of the total spectrum by -1-0.34 % to account 
for the lower sensitivity during observation 1452. In this way, 
the data for both stars are all corrected to a similar PHA-level. 

The spectra were extracted from a region with half- width h 
in the cross-dispersion given hy h - max(l .02, 0.027|/l|) where 
h is in arcsec and A is the wavelength in A. The background was 
extracted from the two rectangular regions between 10"- 40" 
above and below the source spectrum, in the cross-dispersion 
direction. 

The spectra of both stars were fitted with a spline contin- 
uum with knots with a spacing of 5 A; the values of the spline 
at each knot and the associated uncertainties were determined. 
This model was then folded through the response matrix, and 
higher spectral order contamination, although very small, was 
taken into account (see Sect. IA.3l for an estimate of the associ- 
ated uncertainty). From these fits to the count spectra, the ratio 
between both photon spectra was calculated. Thanks to the high 
spectral resolution, this ratio is independent on any errors in the 
effective area. The ratio is shown in Fig.[T] 

2.2. White dwarf model spectra 

The white dwarf models that have been successfully used in 
modelling the spectra of Sirius B and HZ 43A can be classified 
as follows: 

1 . Pure hydrogen models 

2. Homogeneous H-He models 

3. Stratified models with a hydrogen layer atop of a helium 
atmosphere 

We have calculated a set of spectra for these models with 
a range of parameters that are appropriate for Sirius B and 
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Fig. 1. Ratio of the spectrum of Sirius B to the spectrum of 
HZ 43A. data points with error bars: ratio as observed by the 
LETGS. The solid line shows our best-fit model 2 (with a long 
wavelength cut-off of the Lyman pseudo-continuum), and the 
dotted line our best-fit model 1 (short cut-off), as discussed in 
Sect, m The lower panel shows the fit residuals of the observed 
ratio with respect to model 2; the dotted line in that panel shows 
on the same scale the differences between model 2 and model 
1 . Note that the first data point at 50 A is off-scale in this lower 
panel. 



HZ 43A, respectively. The NLTE model atmosphere calcu- 
lations were done wit h the latest version (v. 202) of Tlusty 
( Hubenv & Land 1995 ). The plane-parallel model atmospheres 
were discretised with 100 depth points which span a wide 
range in Rosseland optical depth from 10"^ to 10"^. The whole 
spectrum is sampled with about 7 000 frequencies between 
4 X 10'^ and lO'^ Hz. Hydro gen was represented with a 16- 
level model atom supplemented with a higher superlevel com- 
bining all higher levels up to n = 80. Treatment of level disso- 
lution, pseudo -continuum, superlevel, and line opacity follows 
Hubenv et al ] fl994). Model atoms of He i and He ii incorpo- 
rate 24 and 20 levels, up to n = 8 and 20, respec tively; a descrip- 
tion of the helium atomic data can be found in lLanz & Hubenv 



(I2003h . 
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In most cases the differences between the models in the UV 
and optical part of the spectrum are small, which is no sur- 
prise as the range of parameters has been adopted from previ- 
ous studies that estimated these parameters by fitting the optical 
and UV spectra. However, as the X-ray spectrum can be con- 
sidered more or less as an exponential tail to the full spectrum, 
a tiny shift in effective temperature or surface gravity can have 
a large effect to the X-ray spectrum. For this reason the X-ray 
spectra need to be calculated with care. 

2.2.1. Compton scattering 

A potential important effect for the radiative transfer of X-ray 
photons in the white dwarf at mosphere is C ompton scatter- 
ing. This was first indicated bv ^ Madei ( 1998 ) for the case of 
HZ 43A. Hard X-ray photons scatter many times before they 
reach the surface of the star, and the combined effect of all 
the collisions with relatively cool electrons can cause a sig- 
nificant energy loss for these photons, and therefore a soften- 
ing of the high-energy tail of the spectrum. This effect is most 
noticeab le at wavelengt hs smaller than 100 A. However, re- 
cently Suleimanov et al.l ( l2006i) showed using a more sophisti- 
cated calculation that for both HZ 43A and Sirius B Compton 
scattering can be completely ignored. Therefore, in the present 
calculations we will also ignore Compton scattering. 

2.2.2. Hydrogen Lyman pseudo-continuum opacity 

When we compared our ca l culatio ns in detail with the calcu- 
lations of Beuermann et al. I (l2006l) (see Sect. 5.5! for more de- 
tails), we got a good agreement for HZ 43A but large differ- 
ences for Sirius B in the EUV/X-ray band. A detailed investi- 
gation of the problem showed that the differences are due to the 
treatment of the Lyman pseudo-continuum. 

What basically happens is the following. In the high density 
atmosphere of Sirius B, the density is high enough that a frac- 
tion of the atoms have neighbours so nearby that they are sig- 
nificantly perturbed. The highest energy levels for these atoms 
become partly dissolved, that is an atom ending in them can 
be viewed as ionised, and these le vels are v iewed as partly dis- 
solved, and partly truly occupied. iHummer & Mihalas (.1988 ) 
have considered this phenomenon in detail, and calculated 
well-defined occupation probabilities for hydrogenic levels. 

Consequently, in addition to a traditional continuum that 
extends from the photoionisation threshold shortward and that 
corresponds to a true photoionisation, there is also a "pseudo- 
continuum" that extends longward of the threshold, and that 
corresponds to transitions from a bound lower level to the 
dissolved parts of higher levels. The basic physical process 
is well-defined, because a transition to the dissolved part of 
the level leaves an atom indeed in an unbound state and 
thus the process is a sort of photoionisation, but the essen- 
tial problem is how to formulate the appropriate cross-section 
for the pseudo-continuum. According to t he standard approach 



called "dissolved fraction" that is given through a dissolution 
probability of a fictitious level that would correspond to the 
current freque ncy v (for d etails, refer to the above cited pa- 
pers). Hubenv e t al.l ( 119941) outlined a proof of this assertion, 
but stressed that the proof only applies "close to the photoion- 
isation limit". It is hard to specify exactly how far from the 
edge is the formalism valid, but it was certainly never meant to 
extend hundreds or even thousands of Angstroms away from 
the edge. Therefore, most researchers begun to use ad hoc cut- 
off frequencies for the pseudo-continuum cross-sections, with 
the belief that their actual values do not significantly influ- 
ence modelling results. Some researchers have even introduced 
smoothed cut-off regions (Bergeron, private communication), 
but in any case it should be clearly understood that any cut-off 
is a completely ad hoc concept. 

In fact, a more rigorous treatment of this problem should 
be developed, perhaps analogously to a recent treatment of 
the pseudo-con tinuum cross-sections for neutral perturbers by 
Kowalski & Saumon (2006) and Kowalski (2006), but this was 
not yet done. For the moment, we are left with a necessity of 
using artificial cut- offs. The TMAP co de (IWerner et al.ll2003h 
employed by Beuermann et al.l (l2006h uses a long cut-off at 
243 1 A, while the default for Tlusty, employed in the present 
work, is a short cut-off at 925 A, close to the Lyman edge. 
Lacking a better physical model for the pseudo-continuum, it 
is hard to decide which approach is better. Therefore we use in 
this paper two different sets of calculations: one with a short 
cut-off (referred to here as model 1), the other with a long cut- 
off (model 2). 



2.3. Scaling of the white dwarf spectra 

The model spectra calculated with our code gives the 
Eddington flux Hy (usually expressed in units equiva- 
lent to Wm^^Hz '). The photon spectrum A^^ (in pho- 
') seen at Earth is then given by 



hA 



(1) 



(iDappen et alJI 19871 : iHubeny et al.lll994h . the cross-section is 



formulated through an extrapolation of the traditional pho- 
toionisation cross-section shortward of the edge, and the so- 



where T{A) is the transmission of the ISM and fd = jS 
with R the radius of the star and d the distance to the star. 

For Sirius B, the distance is known accurately with a pre- 
cision of 0.4 % (the HIPPARCOS-based pai-aflax is 0".3792 ± 
0".0016. lHolberg et a l. 1998). Its radius R can be determined in 
principle from the surface gravity g - GM^ /R^, since the white 
dwarf mass Mb is known with an accuracy of 1.5 %. However, 
the typical uncertainty in g is 15 % based on fits to the optical 
and ( extreme) ultraviolet spectra (for example iHolberg et al. 
1998b . Using the currently best available value for the gravi- 
tational redshift (see later), the uncertainty in g is still 12 %. 
Therefore, using this scaling in the form fd - GM^jgcf gives 
a flux uncertainty of at least 12 %. 

There is a simple way out, however We define /' = fdg 
and substitute fd - fig in ([T]i. For a given spectral model, g 
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is prescribed as one of the input parameters of the model and 
therefore exactly known. It is easy to show that 



/' 



(2) 



with flp — 2".490 + 0".004 the photocentric semi-major 
axis dCatewopd & Gatewood' 'l978"), P = 50.09 + 0.06 yr 
the orbital period (van den Bos 1960), and ji^ - Mb/(Ma + 
Mb) = 0.3295 + .0010 the relative mass of Sirius B 



(IGatewood & GatewoodI 119781) . Substituting these numbers, 
we obtain a value off = (2.084 ± 0.019) x 10"'"* ms^^, im- 
plying that we can determine - given the surface gravity cor- 
responding to the model being considered - the flux with an 
accuracy of 0.9 %. 

For HZ 43B, the distance is less accurately known, and we 
use here the observed flux in the optical band to scale the spec- 
trum. The most accurate and best calibrated flux is obtained 
inJieV band: the V magnitude of HZ 43 i s 12.909 ± 0.0017 
dBohhnllioOOl) . iHolberg & BergeronI (l200d ) give the relation 
between the monochromatic flux at 5423 A and the V mag- 
nitude as Fa = 3.804 x 10"'' x 10-°^^^, within about 0.5 %. 
This corresponds to a photon flux of 71.25 photm"^ s"' A"' at 
5423 A. We use this flux to scale all our model calculations, 
and assign a nominal uncertainty of 1 % to it. Note that a sim- 
ilar procedure is harder to apply to Sirius B, as contamination 
by Sirius A enhances photometric uncertainties (see Sect. |2.4t . 



2.4. Constraints on the spectral models 
2.4.1. Gravitational redshift 

The gravitational redshift of Sirius B is now known much 
more accurat e ly usi ng STIS data (vg - 80 ± 5 kms"'. 



Barstow et al.l |2005[) then prev iouslv (for e xample 



1 6 km s~ ' , iGreenstein et all 1 97 ll: 85 + 1 5 km s ~ ' . iHebrard et al 
1999h . This more accurate gravitational redshift, combined 



with the mass derived from astrometric methods allows an in- 
dependent estimate of the surface gravity g = vlc^ IGMb or 



(3) 



using the same fundamental parameters as employed for (|2]i. 
This leads to \ogg (ms^^) = 6.62 + 0.05. In our spectral mod- 
elling, we have discarded any solution where \ogg deviates 
more than 2cr fr om the above value. 

For HZ 43A.lReidl ( l 19961) gives a value ofvg = 30.1 kms 
Kruk et al. ( 2OO2I) estimate the accuracy of this value to be 10- 
15 km s~' . Therefore we limit our models to those cases where 
the gravitational redshift of HZ 43 A is within the 10-50 km s"' 
range. We do this as follow s. The parallax of HZ 43 A is 



known to be 15.3 + 2.9 mas (van Altena et al. 1995). We use 



this to derive the minimum and maximum allowed distance d. 
For a given surface gravity, this distance range corresponds to 
an allowed range for Ug = gd ^[fd|c when we use the accu- 
rately known value of - R^/d^ derived from the optical flux 
(Sect. l2.3T l. We then add this uncertainty in quadrature to a nom- 
inal uncertainty of 15 km s ' in the observed redshift. Solutions 



off'by more than 2cr are discarded. Also whenever the predicted 
range for Ug does not overlap with the 10-50 kms ' range, we 
discard the solution. 



2.4.2. Optical and UV flux of Sirius B 

For HZ 43A, we already scale our spectra to agree with the 
optical flux of this star. Our model spectra for Sirius B must 
also be in agreement with optical and UV flux measurements 
(Table We only consider here monochromatic fluxes, and 
avoid the use of magnitudes as this involves an extra compli- 
cation, namely the convolution with filter transmissions. Care 
should be taken in assessing the uncertainties in those flux mea- 
surements. 

Hebrard et al.1 d 19991) analysed the HST-GHRS Echelle-A 
spectrum of Sirius B. From their Fig. 2 we measure the con- 
tinuum flux at 1302 A as 1.270 x lO"'^ Wm^^s-i A"' or 



(8. 32+ 0.07) X 10"* photons m"^ s"' A"'. Although the statistical 
errors are small, systematic effects are larger as shown below. 



iMack et al 



(119971) have calibrated the GHRS Echelle-A 
spectrometer by comparing spectra with lUE spectra of the 
standard star jj Col. From their Fig. 4c we find that the relative 
calibration accuracy of this instrument between 1250-1350 A 
is 0.84+0.36 percent (the residual r.m.s. scatter between the 
GHRS and lUE spectra). Thus, the flux depends critically upon 
the calibration of lUE. 



Bohlin et al.1 (Il990l) used more than 2700 individual lUE 



SWP and LWR spectra to define the absolute flux distributions 
of the 37 HST standard stars in the wavelength range 1150- 
3300 A. They conclude that the systematic external errors in 
the fluxes are less than 15 percent, while comparison with ANS 
flux measurements demonstrates an internal consistency of the 
lUE spectrophotometry of 2 percent. The basis for the absolute 
flux scale in the UV is given by the spectrum of 77 UMa, which 
h as an uncertainty of 10 % in its absolute flux scale according 



to 



BohUn et all (Il990l) . lUE fluxes may be too low by ~ 10 ' 



In summary, we combine the following eiTors: systematic 
uncertainty lUE scale: 10 %; internal uncertainty lUE specrum 
yu Col: 2 %; relative error GHRS/IUE: 0.84 %; statistical error 
GHRS spectrum Sirius B: 0.8 %. The resulting total uncertainty 
is 10.3 %, obviously dominated by the uncertainty in the abso- 
lute UV flux of T] UMa. We therefore adopt a 1302 A flux of 
8.32 ± 0.86 X 10"* photons m"- s"' A"' . For this wavelength, we 
can neglect any ISM continuum extinction. 



Barstow et al 



( l2005h have analysed STIS spectra of 
Sirius B. From their Fig. 5, we estimate the flux at 3500 and 
4600 A (using the G430L grating), and from their Fig. 4, we 
estimate the flux at 6400 and 6700 A (using the G750M grat- 
ing). The typical statistical uncertainty combined with the un- 
certainty on the STIS flux scale is abo ut 1 %, and we adopt that 
as the nominal statistical uncertaintv. lBarstow et al. 1 12005) ar- 
gue that although STIS has been calibrated to a nominal preci- 
sion of ~ 1 %, the necessary usage of the narrow 52"x0.2" slit 
for Sirius B gives an add itional 4.5 % eiTor, as estimated first 
by Bohlin & Hartigl (Il99 8). This eiTor, however, always leads 
to an underestimate of the flux. Therefore, we add a one-sided, 
positive systematic uncertainty of 4.5 % to the fluxes. 
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Table 2. Measured monochromatic fluxes of Sirius B. 



Wavelength 


flux 


Stat, error 


syst. error 


ref. 


(A) 


(phm-2 s-i A-') 


(%) 


(%) 




250 


19.5 


10 


±10 


a 


700 


20 




±100 


b 


1302 


8.32 X lO'' 


0.8 


±10.2 


c 


3500 


9.41 X lO'' 


1.0 


+4.5 


d 


4600 


6.62 X 10^ 


1.0 


+4.5 


d 


6400 


2.64 X 10^ 


1.0 


+4.5 


e 


6700 


2.25 X 10^ 


1.0 


+4.5 


e 



remaining residu als of the order of 5 %. Also the analysis of 



a EUVE MW,'Holberg et al.7l998!) 
b EUVE LW, Craig et al. ( 1997) 

c HST-GHRS. Hebrar d et al.l (Il99^ 

d HST-STIS G430L. ,Barstow et alj j2005l) 
e HST-STIS G750M jBarstow et al.1 ( l2005h 



In practice, for the five UV and optical flux values we add 
the statistical and systematic uncertainties in quadrature, and 
calculate for each point the contribution to a formal when 
comparing to a set of model fluxes. In this, the one-sidedness 
of the systematic uncertainties in STIS fluxes is taken into ac- 
count. Whenever this (with 5 degrees of freedom) exceeds 
the 2cr upper limit (eitx^ =5+2 VlO) we discard the solution. 



Sing et alj (120021) of a sample of 7 DA white dwarfs shows sys- 
tematic residuals up to 10 %, with a typical scatter of 5 %. We 
therefore assume that the absolute fluxes measured with EUVE 
have a systematic uncertainty of 10 %. 

We constrain our models using the EUVE flux measure- 
ment at 250 A (1.55 x 10""' Wm^^ A"'), taken from Fig. 1 
of iHolberget al.1 (ilQQSh . with a 10 % statistical error mar- 
gin added in quadrature to the systematic uncertainty, to ac- 
count for both the statistical and systematic uncertainties in the 
EUVE spectrum. When the EUVE flux point at 250 A exceeds 
the Icr bound, we discard the solution. 

Finally, we list in Table |2] a constraint from the longest 
wavelength EUVE LW spectrometer, although the error on this 
estimate is large. The main reason for inclusion is that the spec- 
trum at this wavelength is rather sensitive to the foreground ab- 
sorption model, and we want our model not to be in conflict 
with this EUV constraint. 

2.4.4. Interstellar absorption 



Table 3. Adopted absorption parameters for Sirius B and 
HZ 43A. Numbers in brackets are derived from the basic num- 
bers given in the other entries of this table. 



2.4.3. EUVE flux of Sirius B 

The EUV flux of Sirius B is a very s ensitive indicator fo r 



the spectral parameters of this source (Holb erg et al. 



Unfortunately, there is little information publicly available 
about the effective area calibration of EUVE. 

The following information is given in the EUVE guest ob- 
server handbooly. It explains at page 2-15 the calibration pro- 
cedure, paraphrased below. "Because of the lack of standard 
stars in the EUV, the effective area of the spectrometers was 
determined from model spectra for continuum sources. White 
dwarf stars provided the baseline measurements. After spectra 
were extracted, white dwarf model spectra for the object were 
input to a spectrometer simulation program, and the model pa- 
rameters adjusted to produce simulated data that match the ob- 
served spectra as nearly as possible. In most cases, T^ff and g 
were well constrained by other optical and UV spectra. The re- 
sulting models were then compared whenever possible to pho- 
tometric and spectroscopic EUV/soft X-ray observations of the 
same target made with the ROS AT WFC, EXOSAT, HUT, lUE, 
and various rocket experiments. The effective area function was 
then readjusted to bring the measurements into better agree- 
ment." 

However, given the uncertainty of white dwarf fluxes mea- 
sured with previous instruments in both the UV and soft X- 
ray band, as well as the strong dependence of the soft X-ray 
flux on slight changes in the spectral parameters, we believe 
there is significant systematic uncertainty involved; moreover, 
the ground calibration was accur ate to about 2 5 %. A fit to 



the EUVE spectrum of HZ 43 bv lBarstow et al.l (119951) shows 



Ion 


Sirius B 


HZ43A 


Hi column (10^^ m"^) 


0.58 ±0.11 


0.85 ±0.06 


He I column (10^^ m"^) 


(0.045±0.010) 


0.057±0.002 


Hen column (10^^ m"^) 


(0.017±0.004) 


0.039±0.012 


Hi/Hei 


12.8±1.4 


(14.9±1.2) 


Heii/(Hei+Heii) 


0.27±0.04 


(0.41±0.08) 



a EUVE M W.lHolberg et al.l (Il998h 

b EUVE LW. ICraig e t al. (1997|) 

c HST-GHRS.lHebrard e t al.i (ll^9^ 

d HST-STIS G430L. Ba rstow et alj_(l2005h 

e HST-STIS G750M. lBarstow et alj ( l2005l) 



The adopted values of the interstellar absorption towards 
Sirius B and HZ 43A are shown in Table [3] The hydrogen 
column towards Sirius B is the weig hted average of (0.65 + 
0. 16) X 10^2 m-2 (|HebrardetalJll999l from Lya with the HST 
GHR S detector) and (0.52 + 0.14) x 10^2 m ^ dHolberg et al. 



available at http://archive.stsci.edu/euve/handbook/handbook.htmll 



19981 from Lya with lUE detector). Unfortunately Sirius B 
has a low EUV flux such that only upper limits to the He 
column densities are known : < 0.03 x 10^^ m ^ for He i and 
< 0.14 X 10^2 m-2 for Heii dWolff et al.lll999l) . Therefore, we 
adopt the average H i/He i ratio of 12.8 + 1 .4 derived for the lo- 
cal enviro nment of the S t in (|Slavin & Frischi.2007.) . Similarly, 
we follow iHolberg et al.l d 19981) who use the averag e He ii/He 
ratio of 0.27 + 0.04 derived by [Barstow et all (Il997[) . Note that 
the uncertainty in the hydrogen column of Sirius B leads to an 
uncertainty of 34 % in the predicted flux at 700 A. At short- 
her wavelengths the corresponding flux uncertainty is smaller, 
for example only 2 % at the peak of the soft X-ray spectrum 
(150 A). 
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The neutral hydrogen colu mn towards HZ43A (Table [3]l 
has been discussed in detail by iKruk et alj (l2002h . and we use 
their value. F or He i and He ii we use the weighted average of 



Dupuis et al.k l995): Wolff et al. (1999); Barstowe^t al, (1997), 



based on EUVE measurements. 

Due to the low column densities, other ions than those from 
hydrogen and helium can be ignored in the Chandra band: 
we estimate that between 1-180 A metals contribute less than 
0.3 % to the continuum opacity. For our opacity calculations, 
we use the model by Rumph et al. ( 1994). The uncertainty in 
the measured column densities hardly affects our estimated flux 
of HZ 43 A in the Chandra band: at 180 A where the ISM opac- 
ity is highest, the uncertainty in the column densities of H i. He i 
and Hen corresponds to a transmission uncertainty of 0.3, 0.2 
and 1.0 %, respectively. The combined r.m.s. error is 1.1 %. 
Hence, the uncertainty in our model is mainly determined by 
the intrinsic parameters of HZ 43A itself. 

In our spectral modelling, for both stars we allowed the col- 
umn densities for HZ 43 A and Sirius B (only H i) or the ratios 
Hi/He I and Heii/He (for Sirius B) to vary between the +2cr 
limits as given in Table [3] 

2.5. Parameter estimation 

Our model depends on twelve parameters, six for each star, 
namely the effective temperature, surface gravity, helium abun- 
dance or hydrogen layer thickness, and interstellar column den- 
sities of Hi, He I and He ii. Often there are strong correlations 
between these parameters, and we have a number of constraints 
to be obeyed (Sect. lZ4] i. In summary, these constraints are: 

1 . 250 A EUVE flux of Sirius B 

2. Optical and UV fluxes of Sirius B at 1302, 4600, 6400 and 
6700 A 

3. log^ for Sirius B within the 2cr range of the value derived 
from the gravitational redshift 

4. Gravitational redshift of HZ 43A 

5. derived column densities (or ratios) towards both stars 
within their 2cr error limits. 



We define as usual a quantity x given by 



(4) 



where (9, and AO, is the observed LETGS ratio and the associ- 
ated uncertainty of the Sirius B to the HZ 43A spectrum, and 
Mi is the predicted ratio based on our model. We use the data 
points between 50-175 A with a spacing of 5 A as derived in 
Sect. 12.11 Whenever any of the constraints of Sect. 12.41 is vi- 
olated, we formally add io a large number (1000) in order 
to discard that solution. However, as it is more likely that our 
constraints are near the expected value than at their extremes, 
we add for each of the above five constraints a nominal h.x^ to 
(lUl corresponding to the number of standard deviations for that 
constraint. 

We find the best solution using a Monte Carlo method. 
Starting with a broad range of allowed parameters, we draw 



random sets of parameters within that allowed range, and eval- 
uate for each set. Solutions with larger than a thresh- 
old are discarded. After having obtained a sufficient number 
of solutions, we slowly decrease the x^ threshold and simul- 
taneously shrink the allowed parameter space, encompassing 
with some margin all solutions that up to then have been ac- 
ceptable. All acceptable solutions are stored, and after having 
reached the best solution with - ■ we find the errors on 

^ ^ min 

the parameters by finding for each parameter the minimum and 
maximum value for which < + I. We also store each 
acceptable spectrum, so we can also determine for each wave- 
length the range of allowed flux values. 

3. Models to be used for Sirius B and HZ 43A 

3.1. Sirius B 

For Sirius B we have used homogeneous models, which include 
a pure hydrogen atmosphere as limiting case. 

We have also calculated a grid of stratified models for 
Sirius B, but we were not able to obtain successful fits. 
Basically, we constrained the photometric hydrogen column 
to the range of (1.00 - 1.25) x 10~'^ M^, a r ound fl ie value 



of 1.13 X 10 ' ' Mo found bv iHolberg et all (119981) for this 
class of models. The main reason for the failure is that the 
stratified models show a flux deficit of up to a factor of 2-3 
around 50 A as compared to homogeneous models (see also 
Fig. IS; the deficit sets on below 80 A. As this ra nge was at the 



short w avelength end of the EUVE spectrometer. lHolberg et al 



(Il998h were not able to exclude this class of models com- 
pletely. Thanks to the sensitivity of Chandra it is now possible 
to rule out this class of models. 

3.2. HZ 43A 

For HZ 43A we first consider the homogeneous models. 
Barstow et al W995) have put strict upper limits to the amount 
of He in HZ 43 A, based on the limits to the 304 A line of He n 
in the EUVE spectrum. The nominal equivalent width of this 
line derived by Barstow et al. is 0.2±0.1 A, but due to pos- 
sible systematic effects in the EUVE spectrum this cannot be 
regarded as a detection. For their mixed He/H models, they ob- 
tain an upper limit of 3 X 10"^ for the He/H ratio. We have 
calculated a grid of homogeneous models with He/H ratio's be- 
tween and 10"^. Our models with a small ratio such as found 



bv lBarstow et al.l (Il995h yield fluxes in the Chandra band (10- 
180 A) that are 1.4-2.7 % smaller than the fluxes for a pure H 
model, for the same values of Teff and g. It is clear that such 
small differences can be easily accommodated for in a pure H 
model using slightly different values for T^ff and g, which are 
still consistent with the limits from other parts of the spectrum 
to these numbers. We conclude that - at least for our calibration 
purposes - we can safely adopt a pure hydrogen model as far 
as the class of homogeneous models is concerned. 

The other important class of models that include He are 
the stratified models. We have made a grid of models with a 
hydrogen layer mass between lO"''* and 10 M©. All mod- 
els with a hydrogen layer less than 10 M© produce too 
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100 

Wavelength (A) 

Fig. 2. Spectrum of HZ 43 for stratified models with hydrogen 
layers (in M©) as indicated on the plot, divided by the spectrum 
of a pure hydrogen model. The calculation is done for log^ 
(m s"^) = 5 .946 and Tetf = 5055 1 K (the best-fit parameters for 
model 2, see Table |5]l. 



deep He features in t he sp ectrum, consistent with the find- 
ings of Barstow et al. ( 1995b . On the other hand, if the hydro- 
gen layer mass is > 10"" M©, the spectrum cannot be dis- 
tinguished from a pure hydrogen atmosphere in the Chandra 
band. However in the case of an intermediate thickness hydro- 
gen layer of ~ 10"'^ there is a significant eff'ect in the 
Chandra band (Fig. |2|i. For the lowest allowed hydrogen col- 
umn (2 X 10"'^ Mq, see below), there is a ~70 % flux reduction 
around 55 A, diminishing rapidly for larger columns, while at 
shorter wavelengths the flux is higher, peaking for 4x10"'^ Mq 
and then diminishing rapidly in the pure hydrogen limit. 



Table 4. Lyman series of He n in HZ 43. 



Line 


A (A) 


EW (A) 


Lya 


304 


0.14±0.11 


LyyS 


256 


0.05±0.07 


Lyr 


243 


0.09±0.06 



For this reason, we carefully reconsider the constraints to 
the hydrogen layer. From the EUVE spectrum published by 
Barstow et al.l (Il995h (their Fig. 2), we estimated the equiva- 
lent widths of the He ii Lya, /3 and y absorption lines (TableHJi. 
We compared these equivalent widths with our model calcula- 
tions. None of them can be regarded as a detection. Combining 
all three lines, we obtain a best fit hydrogen mass M = 2 - 4 x 
10"'^ Mq. This mass predicts equivalent widths for the He ii 
Lya, j6 and y lines of 0.02, 0.04 and 0.02 A, respectively. The 
99 % confidence lower Hmit is 1.4x10"'^ Mq. While the lower 
limit is rather strict, we cannot fully exclude that the hydrogen 
layer is very thick (which would make the model effectively a 
pure hydrogen model): the pure hydrogen limit is at the 90 % 
confidence upper limit. Another constraint on the thickness of 
the hydrogen layer is obtained from the He n Lyman limit edge. 



Effectively, the edge is slightly shifted due to blending with the 
higher Lyman series lines and occurs near 230 A. We find from 
the EUVE spectrum that the edge is invisible, corresponding to 
an optical depth of less than 2 %. This corresponds to a lower 
limit to the hydrogen mass of 2.2 x 10 Mq. 

We conclude that the thickness of the hydrogen layer in 
HZ 43 is most likely between 2^x10"'-' Mq, with lower values 
excluded but with no solid upper limit to the hydrogen mass. 

Given all this, we consider only stratified models for 
HZ 43A with hydrogen mass > 2 x lO"'-' Mq; as argued be- 
fore, a pure hydrogen atmosphere is a limiting case of this set 
of models. 

4. Results 

Our best fit model 1 (short cut-off) has a of 56.17(52.50), 
our best fit model 2 (long cut-off) has = 49.41 (46.33). The 
numbers in brackets denote the contribution of the Chandra 
data only. With 25 data points and 12 adjustable parameters 
the number of degrees of freedom would be 13, and hence the 
value of is slightly enhanced with respect to purely statis- 
tical noise. However, the actual number of degrees of freedom 
is higher, as several parameters are strongly correlated and the 
best fit is rather insensitive to others (such as the helium abun- 
dance in Sirius B, and in general the interstellar absorption col- 
umn densities). Moreover, our additional constraints also effec- 
tively limit the number of degrees of freedom. Although hard 
to estimate exactly, the true number of degrees of freedom may 
be of the order of 20. 

Fig. [T] shows the observed ratio of the LETGS spectra of 
Sirius B to HZ 43A, together with the best fit models 1 and 2. 
From this figure it is clear that there is some additional system- 
atic scatter in the data points (as the models in the 50 - 170 A 
wavelength range are, as expected, rather smooth). For in- 
stance, the data point at 160 A deviates by -1-3. Oct or -1-3.7 %. 
In this case, some of the systematic effect may be due to the 
fact that this wavelength is close to the edge of the spectrum 
in the -1 spectral order (the physical edge of the detector). For 
other data points, the relative deviations are smaller or less sig- 
nificant. By adding a systematic uncertainty of only 1 or 2 % 
to our ratios, the x^ for the best fit model 2 would reduce from 
49.41 to 29 or 16, respectively, i.e. in the acceptable range given 
the ~20 degrees of freedom. This reduction by a factor of ~2 
in x^ then suggests that we should use Ax^ = 2.0 instead of 
Ax^ -1.0 for the original fits without systematic uncertainties, 
in order to determine the Icr confidence limits on the parame- 
ters. 

We list the best-fit parameters in Table |5] and the spectrum 
at a few selected wavelength in Table|6] The absorbed spectrum 
of HZ 43 is represented with an accuracy of better than 0.5 % 
over the full 43-180 A wavelength range by 



A^,i =exp ( -491.51/^+12.277 

-0.01418/1 -H 11.8 X IQ-^A^) 



(5) 



for model 1 and for model 2 with an accuracy better than 0.7 % 
by 

Na = exp ( -509.12MH- 12.366 
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Table 5. Best-fit parameters of Sirius B and HZ 43A. 



Parameter 


Search range 


Best value" 


Median value* 


Allowed range^ 


Best value 


" Median value* 


Allowed range'' 






Short cut-off Lyman pseudocontinuum 


Long 


cut-off Lyman pseudocontinuum 








(model 1) 






(model 2) 




HZ 43A: 
















Teff (K) 


47000-57000 


51660 


51660 


51420-51880 


51530 


51460 


51240-51690 


\ogg (ms-2) 


5.8-6.3 


6.064 


6.065 


6.031-6.102 


5.893 


5.913 


5.872-5.941 


H-layermass (10"" Mq) 


0.02-OO 


5 


0.4 


0.3-OO 


11 


10 


0.25-OO 


ISM Hi (10-2 


0.73-0.97 


0.85 


0.85 


0.81-0.90 


0.85 


0.84 


0.79-0.91 


ISM He I (10^2 1^-2-) 


0.053-0.061 


0.057 


0.058 


0.055-0.059 


0.057 


0.055 


0.055-0.058 


ISM He II (10^2 m-2) 


0.015-0.063 


0.038 


0.042 


0.029-0.047 


0.038 


0.042 


0.029-0.048 


/, = R^lcf (10-23) 




3.037 


3.037 


3.025-3.050 


3.038 


3.041 


3.031-3.053 


Grav. redshift (km s"') 


10-50 


45 


45 


41-50 


30.0 


31.5 


28.6-33.6 


Sirius B: 
















Teff (K) 


24000-26000 


25360 


25360 


25320-25410 


24970 


24980 


24940-25010 


log^(logms-2) 


6.52-6.72 


6.628 


6.626 


6.624-6.631 


6.624 


6.623 


6.621-6.627 


He/H ratio (xlO-*) 


0-6 





0.2 


0-0.5 





0.6 


0-1.0 


ISM Hi (10^2 m-2) 


0.36-0.80 


0.58 


0.57 


0.50-0.66 


0.59 


0.59 


0.50-0.69 


ISM Hi/He I 


10.0-15.6 


12.8 


13.0 


11.8-14.0 


12.8 


12.4 


11.9-13.9 


ISMHeii/He 


0.19-0.35 


0.27 


0.27 


0.24-0.30 


0.27 


0.25 


0.24-0.30 


/<,=i?2/d2 (10-21) 




4.91 


4.93 


4.87-4.95 


4.95 


4.97 


4.92-4.99 



a Best values (minimum^ ). 

b Corresponding to the parameters of the spectrum closest to the median of all allowed spectra, 
c Corresponding to Ajt'^ = 2 (see text for discussion). 



Table 6. Absorbed fluxes (in photons m ^ s 'A ') of Sirius B 
and HZ 43A at selected wavelengths. 



/1(A) 



Model 1 
Sirius B HZ 43A 



Model 2 
Sirius B HZ 43A 



48 


0.257 


3.98 


0.213 


3.05 


50 


0.416 


5.83 


0.346 


4.54 


60 


2.80 


26.5 


2.39 


21.9 


70 


10.7 


75.5 


9.34 


65.3 


80 


28.3 


160 


25.3 


143 


100 


99.0 


429 


92.1 


400 


120 


192 


772 


184 


738 


140 


253 


1110 


246 


1080 


160 


250 


1400 


243 


1360 


170 


227 


1500 


220 


1470 


250 


19.6 


1780 


18.8 


1760 


700 


15.8 


123 


14.4 


123 


1302 


84000 


2670 


81000 


2670 


3500 


9570 


222 


9540 


222 


4600 


6580 


113 


6590 


113 


5423 


4330 


71.3 


4340 


71.3 


6400 


2720 


44.5 


2730 


44.5 


6700 


2400 


39.0 


2400 


39.1 



6.1 ■ 



, 6.0 



5.9 - 
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Barstow et al. 
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_l I I I L- 



50000 



50500 



51000 51500 52000 

Effective Temperature (K) 



52500 



Fig. 3. Contours of A;^f^ = 2 and h.x^ = 8 in the effective tem- 
perature - gravity plain of HZ 43 A, for model 1 and model 2. 
The cross indicates our best solution for each case. Also shown 
ar e the nominal error ba rs on both parameters from the analysis 
of lBarstow eTaP (l2003h . 



-0.0137Uh-8.8x 10"'^^^). (6) 



Note that (HI)-® should not be used outside this range. 



5. Discussion 

5.1. Parameters of HZ 43 

Our lower limit to the hydrogen mass of 2.5 to 3 x 10"'^ Mq 
is an order of magnitud e higher than the lower limit derived 
by iBarstow et al.l (1 1995b based on the EUVE continuum. We 
could derive this tighter limit because the LETGS covers also 
shorter wavelengths, for which the continuum is very sensitive 
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to the hydrogen thickness (Fig. |2]i. Our best-fit model is in- 
distinguishable from a pure hydrogen model, and even for our 
lower limit hydrogen mass, above 50 A the differences with a 
pure hydrogen model are less than a few percent. 

One of the most recent alternative measure ments of grav- 
ity_and temperature of HZ 43A were given by iBarstow et al 



(only 180 K h igher). For both models, logg is consistent with 
Holberg et al.l(ll99i . 



According to iHolberg et al. 



ity a 

(Eob: 



2003.) based on FUSE observations of the Lyman series lines. 
They obtain values of T^sf = 50 380 + 320 K and log g (m s"^) 
= 5.97 + 0.03. These values are not consistent with our model 1 
or model 2 (Fig. |3]l, but they are close r to model 2. The uncer- 
tainties given by lBarstow et al.l (120031) correspond to the scatter 
between the parameters derived from the individual fits of the 
only three FUSE spectra that were available, hence the nominal 
uncertainty may be quite uncertain by itself. We also note that 
the difference s betw een our best fit model and the model by 
Barstow et al.l ( l2003b are less than 1 % of the continuum level 
in the Lyman series line cores, hence rather sensitive to uncer- 
tainties in scattered light contributions or background subtrac- 
tion. 

Also, for model 1 the allowed range for the gravitational 
redshift is relatively high, given that the best value is 30 + 10 to 
30 + 15 kms"' (see Sect. 12. 4.11 . For the interstellar absorption 
column densities, our model fits do not put strong constraints. 

5.2. Parameters of Sirius B 

Until rec ently, the mos t accur ate parameters of Sirius B were 
given bv lHolbergetaD(ll998 h. These authors used the Lyman 
alpha line obtained by lUE together with the EUVE spectrum 
to constrain the effective temperature and surface gravity. With 
only the lUE observations the effective temperature is known 
within +285 K, but by including the EUVE spectrum this un- 
certainty reduces to ±100 K. 

In a recent paper, iBarstow et al. I (12005') have decreased the 
formal error bars on the effective temperature even further 
down to +37 K by using high-accuracy STIS spectra. However, 
the quote d uncertainty is only the statistical uncertainty, and 
Barstow et al.. (.2005.) argue that the systematic uncertainty on 
these numbers is hard to assess, mainly because there is little 
else to compare with. A major reason of conce rn is the much 



larger effective temperature (25 193 K) found bv lBarstow et al 



(l2005h as compared to the value of 24790 K obtained by 
Holberg et al. 1(11998'). While the surface gravity given in both 
papers is almost equal and consistent within the error bars, this 
temperature difference is > 4cr and causes the EUV flux at 
300 A to increase by a factor of 2.4 according to our own model 
calculations. Although the absolute accuracy of the EUVE cal- 
ibration has its limitations (see Sect. 12. 4.31 1, we believe that a 
factor of 2.4 cannot be easily accommodated for. Moreover, 
Barstow et al . (2005) also show that their best normalisation of 
the G430L spectrum obtained with STIS has systematic uncer- 
tainties larger than desirable. Given this problem w ith the EUV 
flux, we prefer here the older iHolberg et al.l (1 1998b parameters 
with the corresponding error ranges. 

For model 1, the te mperature is clearly h igher (by 570 K) 
than the value found bv IHolberg et al.l (Il998h . while for model 
2 it is consistent with the Holberg et al. value within 2cr 



Sirius B contains a 
small amount of helium so pure hydrogen models are ruled 
out. For homogeneous H/He models, they found tiue/nu - 
(4 ± 1) X 10 Our upper limit to the amount of helium of 
< 1.0 X 10"^ is well below that value. However, the claim of 
detection of helium is based on the non-significant detection 
of possible He ii Lyi5 and Lye lines in the EUVE spectrum, as 
well as on the global fit to the E UVE spectrum . If we evaluate 
our models for the parameters of IHolberg et al.l (1 1998b . we find 
that there should be a deep and sharp He n edge in the model 
near 230 A, with a depth of 13 %. The edge is very broad, and 
reaches half of its maximum depth at 100 A. Clearly, such a 
deep edge is not observed in the EUVE spectrum, and the sys- 
tematic deviations from their best-fit model as shown in their 
Fig. 4 are of the same order of magnitude. In fact, for wave- 
lengths below the He ii edge the EUVE data show even a small 
systematic excess, pointing to a lower helium abundance than 
adopted. We conclude that there is no convincing evidence for 
a substantial amount of helium in Sirius B. 

As for HZ 43 A, our models for Sirius B do not give addi- 
tional useful constraints for the interstellar absorption columns. 



5.3. The cut-off of the Lyman pseudo-continuum 

Based on our fits alone, it is not well possible to distinguish 
between model 1 (short Lyman pseudo-continuum cut-off) or 
model 2 (long cut-off), as both models reproduce well the ob- 
served spectral ratio between Sirius B and HZ 43A (Fig. [T]i, 
albeit with different derived parameters for both stars. When 
we look to those parameters (see previous subsection), it ap- 
pears that the derived temperature for Sirius B and the surface 
gravity of HZ 43A are in reasonable agreement with recent 
literature values only for model 2 (a long cut-off). Model 2 
therefore achieves a better consistency between analyses of the 
UV/optical and soft X-ray ranges. Moreover, model 2 seems to 
match the gravitational redshift of HZ 43 A better 

The most direct test of the Lyman pseudo-continuum is 
provided by the far-UV spectrum of Sirius B. The FUV spec- 
trum was recorded with the far ultraviolet spectroscopic ex- 
plorer (FUSE) on 2002 June 14, using medium resolution 
(MDRS) and the SiC channel covering the range AA916- 
1 100 A. Model 2 with the long cut-off clearly provides the best 
match, though discrepancies up to 10% between the model and 
the FUSE spectrum remain. In particular, the observed spec- 
trum reveals slightly broader high Lyman lines (Ly y and higher 
lines) than those predicted by model 2. The overall continuum 
flux level is however well matched contrary to model 1 . A de- 
tailed analysis of the FUSE spectrum will be presented in a 
separate paper. 

In summary, there seems to be more support for a long cut- 
off (model 2). 
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Fig. 4. Ratio of the spectrum of Sirius B to the spectrum of 
HZ 43A. Data points with error bars: ratio as observed with 
EUVE. The solid line shows our best-fit model 2 (with a long 
wavelength cut-off of the Lyman pseudo-continuum, and the 
dotted line our best-fit model 1 (short cut-off), as discussed in 
Sect, m Note that these fits are based solely upon the Chandra 
LETGS data, not upon the EUVE data. The lower panel shows 
the residuals of the observed ratio with respect to model 2; the 
dotted line in that panel shows on the same scale the differ- 
ences between model 2 and model 1. Note that the first two 
data points below 80 A are off-scale in this lower panel. 
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Fig. 5. Fluxed EUVE spectrum of HZ 43A. Data points: ob- 
served spectrum with EUVE. The solid line shows our best-fit 
model 2 (with a long wavelength cut-off of the Lyman pseudo- 
continuum, and the dotted line our best-fit model 1 (short cut- 
off), as discussed in Sect. |4] The lower panel shows the ratio of 
the observed spectrum with respect to model 2; the dotted line 
in that panel shows on the same scale the differences between 
model 2 and model 1 . 

Table 7. Compari son of model spectra for the same set of pa- 
rameters (those of ^ Beuermann et aL llaooe"). Fluxes are in pho- 
tons m"^ s ' A ' and include interstellar absorption. 



5.4. A comparison with EUVE 

We have obtained fluxed, order-subtracted spectra of both stars 
from the public EUVE archiveH We sampled these fluxed spec- 
tra on a grid with 5 A spacing using a spline fit, and estimated 
the uncertainty on the flux point by looking to the r.m.s. vari- 
ations with respect to this fit in 5 A wide bins centred at the 
grid points (the fluxed spectra from the public archive do not 
contain error estimates). In addition to these statistical uncer- 
tainties, we added systematic uncertainties of about 2, 3.5 and 
5 photons m-2 s"' A"' for A < 150 A, 150 A < /I < 300 A and 
A > 300 A, respectively. These are based upon a comparison 
of the spectra with smoothed spectra on even larger scales of 
~ 50 A. 

The ratio of these fluxed spectra of both stars is shown in 
Fig. |4] It is evident from this figure that the models that we 
found using the observed Chandra LETGS ratios agree very 
well with the observed EUVE data, even beyond the Chandra 
range for A > 170 A. Again, these data cannot help to choose 
between model 1 and 2, although model 2 describes the data in 
the 170 - 250 A range sHghtly better. 

In Fig. |5] we compare the fluxed EUVE spectrum of 
HZ 43A with our models 1 and 2. In the SW band (below 
180 A) the EUVE flux is typically 15 % below our model flux, 
while in the MW band fluctuations up to 10 % occur Note also 
the relatively large systematic fluctuations in both bands of up 
to a few % in flie SW band to 5 % in the LW band. 
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5.5. Comparison witli lBeuermann et a/] \200dli, 



We have 



Beuermann et al 



http://archive.stsci.edu/euve/search.php 



com pared our model calculations with 
by evaluating our model using 
exactly the same parameters as obtained by these authors (their 
Table 2). We show this comparison in Table |7] We have used 
here the long cut-off of the Lyman pseudo-continuum (see 

Sect. 122211. 

It should b e note d that Table 2 and Table 3 of 
Beuermann et al. ( 20061) contain errors, as explained in an er- 
ratum on that paper (in press). This erratum was triggered by 
our present results. Accordingly, we used their updated temper- 
atures of 24897 and 5 1 1 1 1 K for Sirius B and HZ 43 A, as well 
as the updated fluxes (Table 1 of the erratum). 
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Fig. 6. Fluxes of HZ 43 A with respect to the model flux of our 
model 2 (long cut-ofF). Solid line: model 1; circles: results of 
a spline fit to the LETGS spectrum using the old SRON ef- 
fective area calibration based on work in 2000 dash ed his- 
togram: EUVE flux; squares: Beuermann et al. ( 20061) : dash- 
dotted line, triangles: results of a spline fit to the LETGS spec- 
trum using the standard CXC ClAO pipeline. 



For HZ 43A there is an excellent agreement between both 
codes; only at the shortest wavelength listed (48 A), there is a 
small 4 % difference. However, for Sirius B there are large dif- 
ferences. It is striking that at all wavelengths our predicted flux 
is smaller than the flux given by Beuermann et al. (2006), also 
because we used exactly the same interstellar absorption col- 
umn as well as normalisation R^/cfi = 4.877 x 10"^' as these 
authors. We verified that the (unabsorbed) and integrated spec- 
trum of our model obeys with high precision the normalisation 
condition that J F(v)dv - crT^^ with cr the Stefan-Boltzmann 
constant and F{v) the emitted surface flux. As we both use the 
same effective temperature of 24 897 K, the c onclusion must be 
that th e spectrum for Sirius B as calculated bv lBeuermann et al. 
(l2006h is probably not correct. 

Independently, w e compare the ratio of th e spectra of both 
stars as calculated bv ' Beuermann et aL I (l2006l) to the ratio that 
we measured with the LETGS. Again, the measured ratio is 
smaller by on average a factor of 0.941 ± 0.008. The reason for 
this discrepancy is not clear, but we note that our ratio for the 
LETGS spectra is consistent with the ratio obtained from the 
EUVE spectra. 

Remarkably, our own normalisation constant for Sirius B 
for the same v alue of g and T^ff is 7 .3 % higher than the 
value given by iBeuermann et al. I (l2006l) . the main difference 
being that we use a value of R based on the measured gravita- 
tional redshift instead of the spectral modelling derived effec- 
tive gravity. If we would have used our own normalisation in 
Table|71 the discrepancy would have been smaller. 



5.6. Effective area comparison 

As the effective area of the LETGS depends on details such 
as pha (pulseheight) selection, the spectral order, or the width 



of the spectral extraction box, it is not very usefull to give the 
effective area here. Instead, we compare here directly model 
spectra. This comparison is shown in Fig. |6] 

A comparison between model 1 and 2 shows that in par- 
ticular at the shorter wavelengths the differences are large: at 
40 A, model 1 predicts 38 % more flux than model 1. This 
difference is solely due to the adopted value of the Lyman 
pseudo-continuum cut-off that affects the model calculations 
for Sirius B. As we fit both stars together, however, the differ- 
ent spectrum for Sirius B implies then a different solution for 
HZ 43 A, as the ratio of both spectra is constrained by the ob- 
served ratio with Chandra. The difference between both models 
then must be found through a comparison with other data, for 
instance temperatures and effective gravities determined from 
detailed line fitting. 

The differences between the present model 2 and our older 
(2000) effective area based on work by J. Heise and J. Kaastra 
are less than 20 % for ^ > 80 A, and differ by 10 % below 80 A. 
Compared to Model 1, there is a constant offset of about 20 % 
for /I > 80 A; for the shorter wavelengths, the effective area in 
the 2000 version was not based upon HZ 43 but on matching 
blazar spectra from shorter wavelengths, so it is not surprising 
that the differences are larger for shorter wavelengths. 

The model flux given bv lBeuermann et al. I (l2006h IS on av- 
erage 3 % lower than the flux that we derive for model 2, al- 
though the shape agrees within 1 %. 

We also compare our model to the model flux derived di- 
rectly from the observed LETGS spectrum of HZ 43A as pro- 
cessed and modelled using the official CXC software (instead 
of our local SRON software). We used CIAO version 4.0 with 
CALDB version 3.4.2. The spectrum was fitted using a spline, 
similar to what we described in Sect. 12. II and we show in Fig.|6] 
the ratio of this fluxed spectrum to our model 2 with the dash- 
dotted line. Above 90 A, there is good agreement in shape with 
model 1 (apart from a 10 % flux difference), but for shorter 
wavelengths there is a strong dip: between 70 - 90 A, there 
is a relative change of almost 25 % between our models and 
the CXC-based model. The large scale (tens of A) fluctuations 
with an amplitude of a few percent as compared for instan ce to 
our old (2000) calibration and also to IBeuermann et alj (|2006|) 
are not very surprising, as we fudged effectively the effective 
area to get the observed spectrum agree by definition to the 
predicted model. 

Finally, the dotted line shows the comparison of the EUVE 
fluxed spectrum discussed earlier to our model 2. Over the 
wavelength band of the LETGS, the difference with model 2 
is ~ 15 %. 



5.7. Uncertainty of the model flux 

The statistical uncertainty of our model spectrum for both stars 
is shown in Fig. [T] Although not as good as the flux limits in 
the optical band (1 % accuracy), we still can reach an accuracy 
of better than 5 % on the absolute flux of both stars. As we 
will show in another paper (iKaastra et al ] E008bi) . we can even 
reduce this uncertainty to 3-4 % by using additional constraints 
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Fig. 7. Relative uncertainty on our model spectra for Sirius B 
(dashed lines) and HZ 43A (solid lines). All acceptable solu- 
tions with l^x^ < x^^j^ + 2 are bounded by these lines (which 
we scaled to our best solution for each star). Results shown are 
for model 2; the results for model 1 are quite similar. 



at high energies. Here we will carefully asses the possible size 
of systematic uncertainties in our model spectra. 



We follow the same approach as lKaastra et al.l (l2008al) . We 
estimate the known systematic uncertainties and present them 
in Table [8] For each relevant factor, we estimate its magnitude 
for a range of characteristic wavelengths, covering the LETGS 
range between 50-170 A. These numbers are given in the table. 
Then we put each error in one of two specific categories. 

For category "c", the systematic uncertainties at different 
wavelengths are more or less correlated. Example: an error in 
the interstellar absorption column density will lead to corre- 
lated deviations. As a rough approximation, such deviations 
can be approximated by a power law in wavelength. We can 
assess for these "correlated" errors how they will affect the 
normalisation and slope of this power law written in the form 
f{A) = A{A/ 100)" by "fitting" the estimated systematic uncer- 
tainties (all having the same sign) directly to such a power law 
shape. 

For category "u", the systematic uncertainties at different 
wavelengths are uncorrected over the full LETGS wavelength 
scale. A good example are the statistical errors on the Sirius B 
to HZ 43 A ratio. In those cases, we estimated AA from the typi- 
cal slope and normalisation that we would obtain if the system- 
atic uncertainties at different wavelengths had random signs. In 
general, for this type of uncertainty the direct effect on A and a 
is smaller than for correlated errors (category "c"). 

More details on the individual systematic uncertainties are 
given in Appendix A. In summary, the most important factor 
is the uncertainty in the spectral modelling. This can give a 
flux difference at 100 A of 1 % and a difference in slope of 
0.022. However, these values are well within the uncertainty 
range that we found for the spectrum of HZ 43 (Fig. |7]i. 



6. Conclusions 

We have modelled carefully the X-ray spectra of Sirius B and 
HZ 43 A using the most sophisticated spectral models presently 
available. By using the observed spectral ratio of both stars with 
LETGS as constraints, together with some other constraints, 
we have obtained models for each star that predict the flux in 
the 50 - 170 A band to better than 5 %. This allows also to 
calibrate the relative and absolute effective area of the LETGS 
in the same wavelength band with the same accuracy. 

The model spectrum for Sirius B appears to depend sig- 
nificantly upon the assumption made about the cut-off of the 
Lyman pseudo-continuum. This affects our solutions in two 
ways. First, the best fit parameters for both stars depend upon 
this choice, and secondly the absolute flux of both stars at 
the shorter wavelengths depends strongly upon it. Accordingly, 
we have investigated here two classes of models for Sirius B: 
model 1 with a short cut-off of the Lyman pseudo-continuum, 
model 2 with a long cut-off. From our X-ray data alone we can- 
not distinguish between both models. However, it appears that 
only for model 2 (a long cut-off) the derived temperature and 
far UV flux of Sirius B, and the temperature, surface gravity 
and gravitational redshift of HZ 43A agree much better with 
recent analyses of the UV and optical spectra. 

For both stars, we can put tight upper limits to the amount 
of helium or a tight lower limit to the thickness of a hydrogen 
layer In fact, both stars show to agree remarkably well with a 
pure hydrogen atmosphere. The values that we derive for the 
surface gravity and temperature of both stars have an accuracy 
comparable to what has been achieved by detailed line profile 
fitting, apart from the larger differences between model 1 and 
model 2. 
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Interstellar absorption: 

" Uncertainty in the absorption cross sections 

* Uncertainty in the ionisation correction 
' Neglect of metals in the opacity 
Atmosphere modelling: 

Neglect of Compton scattering 

Uncertainties in atmosphere models 
^ Interpolation errors 
Spectral extraction: 

* Background uncertainty 

* Background contamination by source cross-dispersion distribution 
' Higher spectral order subtraction 

' Off-axis vignetting Sirius B 

* Statistical uncertainty spectra 
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Appendix A: Breakdown of systematic 
uncertainties 

A.1. Interstellar absorption 

We have used the interstellar absorption cross sections of 
Rumph et al. ( 19941). Replac i ng the atomic cross sections with 
those of Iverner & Yakovlevl(ll995h gives small differences, as 
listed in column (a) of Tabled 

Also, we use slightly different i onisation fraction s of hy - 
drogen and helium as compared to Beuermann et al.^ ('2006'). 
Replacing our ionisation fractions by those of Beuermann et al. 

(1200 6) yields the differences in colum n (b). 

19941) all 



Further, in the absorption model of iRumph et al 



elements other than hydrogen and helium are ignored. We esti- 
mate this effect by comparing the transmission with and with- 
out metals using the hot absorption model of SPEX with the ap- 
propriate parameters. The differences are very small and shown 
in column (c). 

A. 2. Atmosphere models 

The effects of ignoring Compton scattering has been estimated 
from the study of Suleimanov et al. ( 2006) and are listed in col- 
umn (d). 

Using the adopted model atmosphere for HZ 43A, we es- 
timated the internal accuracy of model atmosphere calcula- 
tions by allowing for different treatment of the hydrogen opac- 
ity (line broadening, extension of Lyman pseudo-continuum) 
and using different discretisation in frequency and depth. The 
cumulative systematic uncertainties are listed in column (e). 
These estimates are similar to the differences found in calculat- 
ing a model with Tlusty and TMAP using the same parameters 
(see for instance Table|7] the column for HZ 43A). 

In our fitting procedure, we calculate atmosphere models 
by interpolating models on a three-dimensional grid (tempera- 
ture, gravity and helium abundance or hydrogen column mass). 
We estimate the uncertainties corresponding to this interpola- 
tion by comparing a few interpolated models with exact calcu- 
lations. We list these in column (f). 



ground subtraction, but the effect can be easily estimated, see 
column (h). 

The higher spectral order contribution to the HZ 43A spec- 
trum is small, less than 3 % for all wavelengths. Adopting con- 
servatively a maximum uncertainty of 30 % on the estimated 
higher order calibration relative to the first order, we deter- 
mined the maximum uncertainty due to higher order subtrac- 
tion as listed in column (i). 

Our observation of Sirius B was taken slightly off-axis with 
respect to HZ 43 A by about 1'. According to the Proposer's 
Observatory Guide of Chandra, the main effect of this off-axis 
angle is due to vignetting of the Chandra mirror At the low 
energies that are relevant for white dwarfs and this off-axis an- 
gle, the vignetting effect is less than 0.5 % and constant as a 
function of wavelength, see column (j). 

Finally, we list the statistical errors on the observed ratio of 
both stars in column (k). 

List of Objects 

'Sirius B' on page! 
'HZ 43A' on page[ 
'PKS 2155-304' on page|2| 
'3C 273' on pagei 
'RXJ 1856.5-3754' onpagei 
'Crab' on page|3 



A. 3. Spectral extraction 

The background towards both stars was determined from boxes 
below and above the spectral extraction region. The back- 
ground of the HRC-S detector is not completely flat. For typical 
HRC-S extraction regions, the background below and above the 
spectrum may differ by 2-3 %. As a conservative estimate we 
therefore can assume that the average background uncertainty 
in the source extraction region is less than 1 %. As Sirius B 
is the faintest of both stars, this uncertainty is the strongest 
for that source. We list the corresponding uncertainty in col- 
umn (g). 

HZ 43A is a very strong source. Inspection and modelling 
of the empirical distribution of photons in the cross-dispersion 
direction shows that at the longest wavelengths the background 
can be enhanced by 10-15 % due to scattered photons of 
HZ 43A. We have not taken this into account in our back- 



